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Abstract 

We discuss constraints on supersymmetric models arising from the relic density measurements of 
WMAP as well as from direct and precision measurements, LEP, b — » S7, (g — 2)„, B s — > fi + u~ . 
We consider mSUGRA models and their extensions with non-universal gaugino masses. We 
find, as commonly known, that the relic density pinpoints towards very specific regions of the 
mSUGRA models: coannihilation, focus and Heavy Higgs annihilation. The allowed regions 
widen significantly when varying the top quark mass. Introducing some non-universality in the 
gaugino masses significantly changes this conclusion as in specific non-universal models the relic 
density upper limit can be easily satisfied. This occurs with models where Mi > M2 at the 
GUT scale when the LSP has a high wino component. Models where M3 < M2 favours the 
Higgs annihilation channel in large regions of parameter space and at large tan (3 also favours 
the annihilation of neutralinos in gauge bosons pairs. We discuss also the potential of direct 
detection experiments to probe supersymmetric models and point out at the main consequences 
for colliders based on the mass spectrum. Our calculation of the relic density of neutralinos is 
based on micrOMEGAs and the SUSY spectrum is generated with SOFTSUSY. 
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Despite the lack of direct experimental evidence for supersymmetry, the minimal su- 
persymmetric standard model (MSSM) remains one of the most attractive extensions of 
the standard model. One of the nice features of the MSSM is that it provides a natural 
cold dark matter candidate, the lightest supersymmetric particle (LSP). Recently, the 
WMAP satellite has measured precisely many cosmological parameters, among them the 
matter density, the baryonic density as well as the relic neutrino density QjJ Ej. From 
these measurements one can infer the relic density of cold dark matter, VLh 2 = .1126j/oi8i 
(at 2a). This experimental measurement severely constrains the parameter space of su- 
persymmetric models j3 IH El- This is particularly true in mSUGRA models where 
most of the analyses have been done so far. Although it has been argued that some 
realisations of mSUGRA are in fact quite generic to more general MSSM models so that 
it is sufficient to restrict to mSUGRA, it rests that most of the interesting scenarios that 
might occur in mSUGRA require a careful adjustment of parameters. For example, the 
Higgsino LSP. Moreover some important scenarios are never realised in mSUGRA like 
for example having a LSP which is dominantly wino. Equality of the gaugino masses at 
the GUT scale (M[ = M' 2 ) leads to M[ = OAM^ at the weak scale thus to a LSP that 
is either mostly bino or a mixed bino/Higgsino state. In fact, in mSUGRA, most of the 
time, at least in what constitutes from the relic density point of view the bulk region, the 
LSP is an almost pure bino. The bulk region corresponds to the low M — M\n region 
of parameter space of mSUGRA models [7[|HJin]. Here M (M 1/ / 2 ) stands for the common 
scalar (gaugino) mass defined at the GUT scale. The recent cosmological data has almost 
ruled out this scenario |U 03 El E3 • This is because the main annihilation channel is into a 
pair of fermions and proceeds through the t-channel exchange of a right-handed slepton, 
the fermions with the largest hypercharge. The annihilation of a pure bino is not very 
efficient since the U(l) coupling, g = e/c, is small and weaker than a SU(2) coupling 
g = e/s. In order to bring down the relic density in the desired range one must appeal to 
specific mechanisms, either annihilation of neutralinos via s-channel Higgs exchange JT] 
or coannihilation of neutralinos with other sfermions ^21- Finally when the LSP has a 
significant Higgsino component, then annihilation into gauge bosons and/or coannihila- 
tion with heavier neutralino / chargino becomes very efficient [Zj , so much so that the relic 
density tends to be below the range of WMAP ^3] • 

Several analyses have shown that at least one of these conditions is satisfied only for 
narrow strips in the parameter space of the mSUGRA model [HI EJ |S] : 

• coannihilation with f\ occur in the low M region and requires m fl »s m x\ (^ ne 
coannihilation region), 

• annihilation via a heavy scalar Higgs into a pair of b-quarks or r leptons is possible 
at large tan (3 (the Higgs funnel region) 

• one finds a neutralino with a high Higgsino content at large values of M when 
the parameter /i becomes small ( the focus point region). This occurs near the 
electroweak symmetry breaking (EWSB) border. 

The strong correlation between the parameters are somewhat responsible for the rather 
contrived constraints coming from the relic density in mSUGRA models. However it was 
also pointed out that when taking into account the uncertainty in the input parameters, 
such as the top quark mass m t , the allowed region widens significantly ^01 IT^j . This is 
particularly true both in the focus point region and in the heavy Higgs annihilation region. 
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In the former, the value of the /i parameter, which drives the Higgsino nature of the LSP, 
is extremely sensitive to m t P^ ll6j. In the latter it is the heavy Higgs masses themselves 
that depend on m t . In the coannihilation region, the top quark mass dependence shows 
up when applying the Higgs mass limit. The direct limit is somewhat relaxed with a large 
top quark mass, therefore we pay special attention to the m t dependence. We confirm 
qualitatively the results of previous analyses [SI IH EJ although we differ in the position of 
the focus point and of the heavy Higgs funnel regions. 

To go beyond the much contrived mSUGRA model while keeping a manageable number 
of free parameters, we introduce non-universality in the GUT-scale relation for gauginos. 
Such non-universality can be found in a variety of models that go beyond mSUGRA. For 
example, SUGRA models with non-minimal kinetic terms superstring models with 
moduli-dominated or a mixture of moduli and dilaton fields and anomaly mediated SUSY 
breaking models [THl H3 1201 El 1221 EH] all feature non-universal masses in the gaugino 
and/or scalar masses. Here we just introduce two free non-universality parameters relating 
the gaugino masses at the high scale, r 12 = Mi/M 2 and r 32 = M 3 /M 2 . 

Basically in non-universal gaugino mass models one finds the same dominant chan- 
nel for neutralino annihilation, however the parameter space for which annihilation into 
leptons, s-channel annihilation near a resonance, coannihilation or Higgsino/wino anni- 
hilation occurs and give a reasonable value for the relic density differs significantly from 
mSUGRA model. Most importantly, we do not have to carefully tune parameters to fall 
in the allowed regions. It is therefore important to explore the relic density as well as 
other constraints from precision measurements on these categories of models in order to 
assess the potential of both future colliders as well as direct /indirect searches for the LSP. 
For one, increasing the ratio Mi/M 2 provides a LSP with a much higher wino content. 
Then gaugino annihilation/coannihilation into gauge bosons becomes the dominant pro- 
cess [21]. Values for the relic density below the WMAP upper bound are found for large 
regions of the M — Mu2 parameter space. On the other hand decreasing the M' 1 /M' 2 
ratio, one can more easily find light neutralinos that can annihilate via a s-channel Z or 
light Higgs pole provided /i is small enough. This was discussed in details in [23 ED 121] • 
A shift in the ratio M 3 /M 2 also induces significant changes in the prediction for the relic 
density. Indeed, the parameter M 3 has a large impact on the spectrum calculation, driv- 
ing the masses of all coloured particles but also changing the relation between jj, and Ma, 
two parameters that are critical for the estimate of the relic density. Thus, decreasing r 32 
introduces noticeable shifts in the region of parameter space where neutralinos annihilate 
via a heavy Higgs pole. In particular this annihilation channel can be relevant at low 
values of tan/3. Furthermore one finds much lighter squarks than in mSUGRA models 
favouring the coannihilation with squarks. 

Although we focus mainly on the relic density constraint, we also take into account 
the experimental constraints set by direct searches at colliders as well as those coming 
from precision measurements, including b — ► sj, (g — 2) M , B s — > fi + fi~. As concerns the 
limit from the muon anomalous magnetic moment, there have been many changes both in 
the theoretical estimates as well as improved experimental bounds in the last years j2H]- 
In particular, the standard model prediction depends on how the lowest-order hadronic 
contribution is extracted from the data, either using e + e~ or r data [29 . We quote both 
the e + e~ estimate or an average of all results. We think it would be overly conservative 
not to include at all this constraint. After all the upper allowed value for 5a ^ has been 
stable over the many refinements of both theoretical and experimental results. 
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When dealing with constraints on the mSUGRA parameters, one crucial step is the 
evaluation of the physical spectrum corresponding to a given mSUGRA model. This is 
based on the renormalization group equations (RGE) that specify the evolution of param- 
eters from the high scale to the low scale as well as on the inclusion of radiative corrections 
to the physical masses of sparticles. This introduces some theoretical uncertainties in the 
prediction of the spectrum and can seriously affect the prediction of observables at low 
energies and in particular of the relic density. Here we will present results obtained with 
S0FTSUSY1.8.6 [HO]- A detailed analysis of the difference between codes for the prediction 
of the relic density as well as a more detailed study of the theoretical uncertainties in 
the evaluation of the spectrum and their implications for the relic density will appear 
elsewhere j2U E2] • 

The relic density measurement is not the only way to probe supersymmetric dark 
matter in space, many experiments are underway and even more have been planned to 
directly search for dark matter. These experiments are effectively sensitive to either the 
spin independent (scalar) or spin dependent neutralino proton scattering cross section. 
The neutralino proton cross section proceeds through Z/Higgs or squark exchange. Thus 
a large cross section requires either light squarks and/or significant coupling with the 
Z/Higgs. In mSUGRA models this means that large cross sections are typically found in 
the focus point region. The current series of experiments have reached a sensitivity on 
the scalar cross section at the level cr^o_ p « lCT 6 pb 34 . The only report of a positive 
signal by DAMA[35] has not been confirmed by others 2 . To really probe mSUGRA models 
however, one still has to wait for improved detectors such as CDMS2 or EdelweissII. CDMS2 
has recently published much improved upper limits [US]- By the time one reaches ton-size 
detectors, such as ZeplinIV[37 j; Genius |HH| and Xenonj^nj a signal could be found for 
a significant fraction of the parameter space of the mSUGRA model. At the same time 
many detectors that are mostly sensitive to the spin dependent part of the neutralino 
nucleon cross section have started taking data, for example Naiad[40 , Simple[41j and 
PICASS0[42j. Upgrades of these detectors are planned and eventually they will reach a 
sensitivity between 10~ 5 — 10~ 4 pb. 

Our calculation of the relic density is based on micrOMEGAsl.3 E3j, a code for 
the calculation of the relic density of dark matter which includes all annihilation and 
coannihilation channels. Although all cross sections are calculated exactly only at tree- 
level, some important loop effects are taken into account: loop corrections in masses as well 
as in vertices, loop corrections for the Higgs widths are also included. The pole masses 
and all mixing matrices are obtained from a code that calculates the supersymmetric 
spectrum. Additional SUSY corrections to the hbb vertices (the Amj corrections ) are 
calculated independently. Although the radiative corrections to neutralinos, charginos 
and sleptons are rather small, they can play a crucial role in the case of coannihilation. 
Other constraints from precision measurements are included as well in the code: b — > sj, 

(g - % , b s -> fi+fi~. 

The paper is organized as follows: we first review the model parameters in section 1 
and the direct and indirect constraints in section 2. In section 3 we present results in the 
mSUGRA scenario taking also into account uncertainties in the top quark mass. We first 
summarize the constraints form colliders and precision measurements before discussing at 
length the relic density constraints. As a first step towards a full analysis of the MSSM we 

2 This discrepancy might be explained using different velocity distributions for light dark matter par- 
ticles insi- 
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then introduce non-universality in the gaugino sector in section 4. The potential for both 
the spin-dependent/spin- independent direct detection experiments for the mSUGRA and 
non-universal mSUGRA models will be discussed in Section 5. Our results are summarized 
in the conclusion. 



1 Model parameters 

The parameters of the MSSM Lagrangian that are necessary for the calculation of the 
various processes are obtained using the renormalization group equations (RGE's) starting 
with the input parameters of mSUGRA models defined at the GUT scale, Mq the scalar 
mass, Mu 2 the common gaugino mass, Aq the common trilinear coupling, tan/3 the ratio 
of the Higgs vev's defined at the weak scale and sgn(fi), the sign of the Higgs mixing 
parameter. The convention used are those of the SUSY Les Houches Accord [45] . we only 
repeat here the convention for the neutralino mass matrix: 



Mi -M z cos (3 sin 6 W M z sin (3 sin 6 W \ 

^ M 2 M z cos (3 cos 6w —M z sin (3 cos 6w 

—Mz cos (3 sin 6 W M z cos (3 cos 6 W —(/, 

\ Mz sin (3 sin 9w —Mzsm(3cos9w ) 

(1) 

Here Mj, the gaugino masses, as well as other parameters entering the matrix are defined 
at the electroweak symmetry breaking scale. The nature of the LSP, Xii which is a linear 
combination of bino B, wino W and the two Higgsino states Hi^, is a crucial parameter 
in the evaluation of the relic density; 

x\ = N n B + N 12 W + Aq 3 #i + N lA H 2 (2) 

where N is the neutralino mixing matrix. By inspection of the neutralino matrix one sees 
that if Mi << M 2 ,fi the lightest neutralino is almost a pure bino, if /i << Mi,M 2 it 
becomes dominantly Higgsino and when M 2 << Mi, [i the neutralino is dominantly wino. 
In mSUGRA models, one always finds that Mi pa 0.4M2 at the electroweak scale, moreover 
in most cases M 2 < /i. Thus the lightest neutralino is mainly a bino. This also leads to 
the usual mass relation between the LSP and the lightest chargino m^o pa .5m~+. Within 
the mSUGRA model, there also exists a possibility for a LSP with a much larger Higgsino 
component, this occurs in the focus point region at large Mq. There, the parameter [i 
decreases very rapidly until eventually /i 2 < and the EWSB cannot take place. When 
\l < Mi,M 2 , the mass of the LSP is driven by \x and one has m^o pa m^o pa The 
coupling of the neutralino to the Z depends on the Higgsino fraction and is therefore 
enhanced in the focus point region. The coupling of the neutralino to the light Higgs 
also requires some Higgsino component and is therefore also enhanced in the focus point 
region of mSUGRA. 

The parameter /i is extracted from the minimization of the potential at the electroweak 
scale. At tree-level, 

_ m? Hl -m^tan 2 /3 1 2 

~ tan 2 f3 - 1 2 Z [6) 
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where m Hl ,m H2 are the two soft scalar masses in the Higgs potential. Thus the value for 
pi strongly depends on the RG evolution of these scalar masses. In mSUGRA one finds 
typically large values for \x as compared to the gaugino masses. In more general models 
defined at the GUT scale one introduces new parameters at the GUT scale which can 
change the relation among parameters at the weak scale, it is then possible to obtain a 
LSP with a significant wino or Higgsino component. 

In models with non-universal gaugino masses we introduce two extra parameters that 
characterize the amount of non-universality, 



« J 1GOT ^3 GUT /.s 

ri2 = r 32 = — (4) 



Mi 


\gut 


M 2 


\gut 



M 3 


\gut 


M 2 


\gut 



where M^gut = M[ are the gaugino masses defined at the GUT scale. We always define 
M' 2 = M 1/2 . 

Relaxing the gaugino masses universality, one can have models where the LSP can 
have a large wino component. For example, ryi > 1, leads to a LSP that has a significant 
wino component. LEP data then implies that m^o > which in turn favours large 
annihilation cross sections for neutralinos into gauge bosons. This also implies that the 
chargino can be almost degenerate with the LSP, m^o ps hence the important 

a-1 K,\ 

contribution of coannihilation channels. 

The parameter r 32 has a strong influence on the weak scale parameters. Not only does 
the gaugino mass M 3 determine the coloured sector but it also enters indirectly the RGE 
equations for Mu x , Mh 2 • In particular, this implies that both \i and Ma also depend on 
M 3 and can be significantly shifted as compared to the unified case. Taking all scalar 
masses equal at the GUT scale and with tan/5 = 10, A = 0, we obtain an approximate 
solution for /i 2 at the weak scale, 

fi 2 ~ -0.5M| + c M 2 + 2.6M' 3 - 0.2M' 2 , + 0.005M'i + 0.034MgM( 

+0.22M£M£ + 0.007A^M( (5) 

The coefficient c depends critically on the top mass at the SUSY scale, for m t = 175 GeV, 
Co ~ 0.11. Co can get very small and even negative which defines the focus point regions 
where \x becomes very small. Clearly then because of the partial cancellation between 
the term in M3 and the one in M' 2) when r 32 < 1, \x is smaller than in the unified case. 
This once again leads to a LSP with a higher Higgsino component and more degenerate 
with the lightest chargino. The former implies more efficient annihilation, the latter more 
efficient coannihilation. 

The pseudoscalar mass also depends strongly on M 3 . An approximate solution to the 
RGE equations under the same condition as above leads to 

M\ sa -M\ + 1.08M 2 + 2.6M'3 + 0.28M'J + 0.22M^M^ + 0.03M' 2 . (6) 

When r 32 < 1, the pseudoscalar mass is lowered as compared to the universal case and one 
can find a pole in the neutralino annihilation channel even for an intermediate value of 
tan/3. For example, when r 32 = 1/2 , m| = .4M| and M ps .7M 2 , we find Ma ~ 2m^o. 
In the universal case such a solution can be found only at large values of tan (3. This will 
open up new regions of parameter space where neutralino annihilation proceeds through 
an exchange of a Higgs in the s-channel. 
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Introducing non-universal Higgs masses would in many ways reproduce the features of 
the non-universal gaugino masses scenarios in the sense that both fi and Ma are directly 
related to the Higgs scalar masses at the EWSB scale. For example, lowering M 3 is 
equivalent to taking Mh 1 larger than other scalar masses. Then it is possible to obtain 
values for \x and Ma that are lower than in the universal case leading to more efficient 
annihilation/coannihilation of neutralinos. We do not discuss these scenarios any further, 
they have been discussed in |46j. 

In this paper we do not use the approximate solutions Eq. which only serve as 
a guide, but obtain the electroweak scale parameters directly by solving completely the 
RGE's. This is done by a call to SOFTSUSY, a program that solves the RGE equations 
and also calculates the radiatively corrected superparticle masses and couplings. These 
higher-order corrections are essential to precisely evaluate the relic density both in the 
coannihilation regions where the exact mass splitting is a crucial parameter and near a 
s-channel resonance. 

2 Constraints on supersymmetric models: Direct lim- 
its, Q/i 2 , (g - 2)^, b -> S7, B s -> fj+fi' 

Direct limits from LEP 

For the chargino mass, we use the LEP2 bound, m+ > 103.5GeV except when there 
is near degeneracy with the lightest neutralino (AM = O(l)GeV), then m+ > 91.9GeV 

HE]. 

For the neutralino sector we do not use directly the lower limit on the neutralino 
mass quoted by the LEP experiments, m^o < 59.6GeV since this value implicitly assumes 
universal gaugino masses at the GUT scale, Mi = M 2 jlOj. Rather, we use in addition to 
the constraint on the chargino mass, the processes e + e~ — > X1X2? XiX3- The associated 
production can somewhat improve on the constraints set from the chargino pair production 
in scenarios where Mi < M2 (SHI EE] in particular in the region with small fi. In our scans, 
we implemented the upper limit from the L3 experiment on a(e + e~ — > X1X2 + X1X3 ~^ 
fi~) PP7| . The radiative process e + e~ — > X1X1I an d the invisible width of the Z, 
T Zinv < 3MeV, can in principle also constrain the neutralino sector but these can play a 
role only in models with light neutralinos Even then, present limits are such that 
these processes do not help much in reducing the allowed parameter space. 

For selectrons, a limit of mg R > 99.4 GeV can be set when m^o > 40GeV while for 
other charged leptons the lower limits of > 96.5GeV and m f > 92.5GeV apply [48j . 

In models where the pseudoscalar mass is heavy, the limit on the Higgs mass from 
LEP2, vtlu > 114.4 GeV, applies. However we have imposed the limit m# > 111 GeV 
to allow for theoretical uncertainties. The Higgs mass is calculated by the code for the 
SUSY spectrum calculation, here S0FTSUSY1.8.6, and includes the two-loop corrections 
of Feynhiggs [21] as well as the two-loop yukawa corrections jH2j. In models with a 
light pseudoscalar the above LEP2 constraint is relaxed, the absolute bound run, Ma > 
91.6 GeV is obtained when sin(a - /3) « [E3 ES|. 
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Relic density of neutralinos VLh 2 

The supersymmetric models with a stable neutralino must be consistent with at least 
the upper limit on the amount of cold dark matter. We consider both the 2a limit from 
WMAP, 

0.094 < Qh 2 < 0.129 

as well as only the upper limit since there could be some other contribution to the cold 
dark matter. Our calculations of the relic density is based on micrOMEGAsl.3, a program 
that calculates the relic density in the MSSM including all possible coannihilation channels 
|44| I43j . Within the context of a model defined at the GUT scale such as universal/non- 
universal SUGRA that we consider here, we first need to calculate the MSSM parameters. 
To do this we use a renormalization group evolution code, SOFTSUSY to evaluate the su- 
persymmetric mass spectrum as well as the mixing matrix elements following the conven- 
tion of the SUSY Les Houches Accord ■ These parameters are then used together with 
tan f3 and fi to calculate all cross sections for annihilation of neutralinos or coannihilation. 
Although the mass and mixing parameters include one-loop corrections, this procedure 
is of course not a complete one-loop computation of the relic density. However the most 
important one- loop effects are taken into account. In particular the mass difference be- 
tween the NLSP and the LSP are computed correctly. This mass difference is crucial in 
evaluating the relic density in the coannihilation region. Indeed in computing the effective 
annihilation cross section, coannihilation processes are suppressed by a Boltzmann factor 
oc exp _AM//T f where AM is the mass difference between the NLSP and the LSP and Tf 
the decoupling temperature (Tf « Mlsp/25). Then even a small shift in the masses can 
reflect significantly on the relic density. For the Higgs width the running of the b-quark 
mass including three loop corrections as well as the Ami, corrections are included 

Muon anomalous magnetic moment (g — 2) M 

The latest experimental data on the g—2 measurement using /n~ |28j . brings the average 

to 

a^ xp - = 11659208 ± 6 x 10~ 10 (7) 

The quantity a M includes both electroweak and hadronic contributions and is still 
subject to large theoretical errors, 

theo. _ QED , weak , HAD , LBL , HAD/NLO 

The largest uncertainty in the calculation of arises from the lowest order hadronic con- 
tribution. We take an average of four recent estimates: the value reported by Jegerlehner 
a had = ( 68g9 ± 5g ) x 1Q -ii ^ 56 j that inc i uc ies BES and CMD-2, the value extracted from 

inclusive data by Teubner et al. a^ ad = (6831 ± 59 ± 20) x 10 -11 [57J, the latest estimate 
by Davier et al. based on e + e~data a^ ad = (6963±62 e2 ,. p ±36 4 / l ) x 10~ n |2n] and finally the 
result of Yndurain, a h ^ ad = (6935 ± 50 exp ± 10 ± 5th) x 10- n [58j. The latter two authors 
also report another estimate that is based on r + r~ data. The values extracted from r 
data are consistent with each other but differ from the values obtained from e + e _ [29 . 
For the hadronic light by light contribution we take a^ BL = (+80 ±40) x 10 -11 jSH] which 
is based on the calculation of Ref. jHOl EH E2J an d includes a large theoretical error. 

We get, after averaging the different hadronic contributions from e + e~ data alone, 

5a^ = a^- - af eo - = (33.8 ± 6| cxp . ± 11.5| theo .) x lO" 10 (8) 



7 



while including also the r estimate we get 

5a M = a^ xp ' - a^ eo - = (28.5 ± 6| exp . ± 11.5| theo .) x 1(T 10 (9) 
Adding the errors in quadrature leads to the 2a allowed range, for the e + e~ data alone, 

7.8 < <5<v x 10 10 < 59.7 (10) 
while including r data in the average, the value moves closer to the standard model, 

2.5 < <5cv x 10 10 < 54.4 (11) 

Note that our central value is slightly higher than the one reported by Davier [2H] 
as we have averaged all estimates of the lowest-order hadronic contribution. We will in 
general only show contour levels for 5a M . We remark that in the last couple of years there 
has been many improvements to the computation of the standard contribution to (g — 2)^, 
and the allowed area has shifted significantly. Still the upper bound has been rather stable 
and therefore should be considered as a conservative constraint on the supersymmetric 
model. 

b — > S7 

Our calculation of the branching ratio for b —>■ sj takes into account next-to-leading 
(NLO) order terms, bremstrahlung and some non-perturbative effects |631 IM] loTj. Our 
standard model value (with scale parameters set at m^), gives Br(b — > 57) = 3.72 x 10 -4 
while the scale and other parameter uncertainty (a s , CKM ...) are about 10% [1H|. The 
Wilson coefficients for the charged Higgs contribution are evaluated at NLOjHH] while 
SUSY contributions at leading-order are included as well as the important large tan/3 

effects inn] 

The experimental value for the branching ratio is extracted from an average of the 
Babar[68 CLEO jEH BELLE [ZD] and ALEPH [H] measurements 

Br{b -»• 57) = 3.34 ± .38 x 10~ 4 (12) 

and assumes a fully correlated theoretical error |72j. 

We require that after allowing for the (scale) uncertainty in the theory calculation the 
result must be within 2a of the experimental result, Eq. El Since the theory uncertainty 
is roughly constant over the SUSY parameter space we have allowed for a conservative 
fixed uncertainty of 10% independently of the SUSY parameters. Thus in effect we require 
the theory prediction to fall within the range 

2.25 < Br{b try) x 10 4 < 4.43 (13) 

The CDF experiment at Fermilab has obtained an upper bound on the branching ratio 
Br(B s -»■ fi+fjr) < 9 x 10~ 6 [73] and should be able to reach Br(B s -> < 2 x 10~ 7 

in Runll. In the SM, this branching ratio is expected to be very small (~ 3 x 10 -9 ). In the 
MSSM, SUSY loop contributions due to chargino, sneutrino, stop and Higgs exchange can 
significantly increase this branching ratio. In particular, the amplitude for Higgs mediated 
decays goes as tan /3 3 and orders of magnitude increase above the SM value are expected 
for large tan/3. Our calculation is based on [21] and agrees with [ZE]- The Amj, effects 
relevant for high tan/? are taken into account. 
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3 Results in mSUGRA 



We have performed scans over the mSUGRA parameter space using the code S0FTSUSY1.8.6 
We first present results in the M — M\m plane for different choices of Aq and tan/3 and 
with fixed values of the top quark mass. We discuss constraints from precision measure- 
ments before concentrating on the relic density constraints. We consider only models with 
H > as they are the ones favoured by the (g — 2) M measurement [ZS]. We note from 
the onset that there are theoretical uncertainties in the evaluation of the spectrum of the 
MSSM within mSUGRA models ^Bj- These can have a large impact on the predictions 
for the relic density [3T]- Here we fix mb(rrib) = 4.23 GeV and choose values of m t within 
the la range, m t = 174.3 ± 5.1 GeV [77j. 3 

3.1 Precision measurements and direct limits 

The mass of the lightest Higgs is basically set by the top and the stop quark sectors. 
Heavier masses are found in regions with large tan (3 and large mixing in the stop sector. 
One expects a sensitivity to the input value of m t , as well as to the value of Aq which sets 
the parameter A t that drives the mixing among the stop quarks. Large mass splitting 
between the stops is difficult to find in the low Mq — Mu2 region of parameter space. 
In this region the Higgs tends to be too light as seen from the contours of constant 
Higgs masses = 111, 114 GeV in Fig. ^ An increase in the top quark mass also 
increases the Higgs mass. This has particularly a large impact on models with small 
values of tan j3. For example, for tan (3 = 5(10) and for small M , the allowed region shifts 
from Mi/a > 417(271) GeV for m t = 175 GeV to M 1/2 > 332(230) GeV for m t = 179 
GeV. For tan (3 = 10, the Higgs mass limit alone translates into a lower bound on the 
LSP m^o > 105(87.4) GeV as well as on the lightest chargino m~+ > 194(160) GeV 
respectively. Of course lighter charginos, in fact down to the direct limit set by LEP, are 
allowed for larger M which correspond to heavier sleptons, but these will in general not 
be relevant from the relic density point of view as we will discuss next. A decrease in the 
top quark mass conversely tightens the constraint on M1/2, for example the m t — 172 GeV 
and rrih = 111 GeV contour roughly coincides with the m t = 179 GeV, = 114 GeV 
contour of Fig. EJi. At large tan f3, the Higgs mass increases as well, so that the constraint 
on both Mo and M1/2 is relaxed. 

The shift in the Higgs mass contours as one varies \Aq\ < 2000 GeV is quite significant 
as displayed in Fig. First note that imposing Aq = in mSUGRA, leads to a negative 
value for A t at the electroweak scale. Since the dominant correction to the Higgs mass 
depends on the mixing in the stop quark sector, oc A t + \ij tan/3, larger Higgs masses 
will be found for Aq < where \A t \ is large. Indeed in Fig. ^p, one sees that with 
Aq = — ITeV, the Higgs mass increases loosening the constraints on the low M — M1/2 
corner of parameter space for tan f3 = 10. Already at tan (3 = 35, one finds rrih > 111 GeV 
over the full parameter space. Conversely, constraints are tightened for A = 2 TeV. Note 
that the region where one can find a consistent solution to the RGE's depends on m t , tan f3 
and on the value of Aq. Hence some of the contours stop abruptly in FigGJwhen one meets 
a non-physical region. Furthermore, the region where the LSP is neutral also shifts with 

3 New results from the Tevatron indicate a higher value for the top quark mass m t = 178 ± 4.3 80 . 
The new central value falls within the range we have considered. 
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Figure 1: a) Contours of rrih = 111 GeV (black) and rrih = 114 GeV (pink/grey) in the 
Mq — Mi/2 plane for \i > 0, A = and for tan (3 = 5 (dot), tan/3 = 10 (full), tan/3 = 35 
(dot-dash). Both contours for m t = 175 GeV and m t = 179 GeV are displayed. For 
each set of parameters, the contour for the heavier m t is on the left. For tan/3 = 35, 
only the contour for m t = 175 GeV is displayed. Some mh,m t values are attached to 
the contours for easy reference, b) Effect of the trilinear mixing for tan/3 = 10, Aq = 
(black) A = —ITeV (green /light grey), A = 2TeV (pink) on the rrih — HI GeV contours 
(full) and Br (b — > sj) = 2.25 10~ 4 contours (dash), c) Contours for rrih — 111 GeV (full), 
tan/3 = 35 with A = (black), A = ITeV (pink) and for Br(b -> s 7 ) = 2.25 10" 4 
(dash) with A = (black), A = —ITeV (green/light grey). 
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tan/9, m t and/or Aq mainly because of the mixing in the r sector as will be discussed in 
the next section. 

The rare processes b — > 37, B s — ► as well as 5a M all provide large supersymmetric 

corrections at large values of tan (3 especially when the sfermions are light. The branching 
ratio for b — > 37 depends mostly on the squark and gaugino/higgsino sector as well as 
on the charged Higgs. Even heavy squarks do not completely decouple and one can get 
substantial corrections to the SM branching ratio. The chargino exchange diagram gives a 
negative contribution relative to the SM one (for \x > 0) and the branching ratio for b — > 57 
drops below the allowed band in the low Mq — M1/2 region. Furthermore there is a strong 
A t dependence from the mixing in the stop sector. For tan/3 = 10 and m t = 175 GeV, 
the b — > 57 constraint is not as severe as the constraint from the Higgs mass as long as 
Ao>0, Fig. For Aq = 2 TeV, the b — > S7 constraint is satisfied over the full parameter 
space. For a large negative value for the trilinear coupling, A = — 1 TeV, the b — > 57 
constraint dominates. As one increases tan/3, for example tan/3 = 35, the dependence 
on Aq is even more striking both for the Higgs mass and for the b — > 57 branching ratio, 
especially that the shifts on the allowed parameter space work in opposite directions. For 
Aq = 2 TeV, the Higgs mass constraint is the most severe since the b — > sj plays no 
role while for Aq = — 1 TeV, the b — > 57 constrains values of M1/2 all the way to 450GeV 
whereas the Higgs mass constraint has no impact, see Fig. ^p. For Aq = 0, it is the b — > S7 
constraint that dominates. In the region of light sleptons which will be relevant from the 
relic density point of view, the lower bound on M1/2 = 345(435) GeV for tan/3 = 35 and 
A = 0(-l) TeV implies m^o > 139(181) GeV and m~+ > 261(347)GeV. As one further 
increases the value of tan /3 one retains the same features: only the size of the forbidden 
region at low Mq — M1/2 increases. 

The branching ratio for B s — > is predicted to be much below current experimen- 

tal limits at low tan/3, however it increases rapidly with tan/3. There the chargino/stop 
exchange contribution dominates and it is largest in the low Mq — M1/2 region [23- This 
constraint depends somewhat on the top quark mass which influences the spectrum in the 
stop sector. The branching ratio is also dependent on the value of the mixing in the stop 
sector. Negative values of A , which as we have discussed above decrease \A t \, strengthen 
the constraint. Eventually B s — ► will become an important constraint, this occurs 

typically for larger values of tan/3 that we have considered [ZHj- Among the cases studied 
here it is only for tan /3 = 50 that in Runll at Tevatron one could put restrictions on the 
allowed parameter space, for example for the case Aq = —1000 GeV shown in FigJSl 

As for (g — 2) M , the allowed parameter space depends critically on how much of a 
deviation from the standard model one imposes. A significant contribution to <5a M is 
obtained with light charginos and light sleptons. The upper limit on Sa^ then rules out 
the lower left hand corner of the Mq — Mxii plane. Larger deviations from the standard 
model are predicted as one increases tan /? but this constraint does not play an important 
role once the limit on the Higgs mass is taken into account. On the other hand, requiring a 
non zero deviation from the standard model implies not so heavy sleptons and neutralinos, 
thus constraining the large Mq — M\ji region of parameter space as well . This region 
is more important at small values of tan/3 where, as just mentioned, one tends to get 
predictions consistent with the standard model, as displayed in Figj2] for the contour 
Sa^ = 5. x 10~ 10 . Introducing a negative value for Aq does not have a large impact on 
the prediction for 5a M as it is the combination A^ + \i tan (3 that comes into play in the 
smuon mixing. For intermediate to large tan/3, the mixing is dominated by /itan/3. 
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To summarize the collider constraints in mSUGRA type models, the mass of the Higgs 
is the main constraint in the low M — M1/2 region for values of tan (3 < 35 and is dependent 
both on the value of the top quark mass and the trilinear mixing. For tan (3 > 35, the 
b — > S7 constraint takes over, at least for Aq not too large or for Aq < 0. If one takes into 
account the lower limit from <5a M one ends up with an allowed band in the Mq — M1/2 plane. 
This band moves towards higher values of M1/2 and increases in size with increasing tan /3. 
We stress again that the upper limit on Sa^ is rather conservative but the lower limit is 
still subject to debate. In the following we will display the contours for fixed values of 



3.2 Relic density 

Within the context of mSUGRA models, one needs a special tuning of parameters to be 
in agreement with the very precise measurement on the relic density by WMAP. We first 
remark that the so-called bulk region (at low My 2 and Mq) that used to be one of the 
favoured regions [7J |H] has considerably shrunk. In the bulk region the LSP is almost 
a pure bino. As such, the LSP couples mostly to fermion/sfermions with the largest 
hypercharge, that is the right-handed sleptons. The dominant annihilation channel for 
neutralino annihilation is into a pair of fermions via the t-channel exchange of a right- 
handed sfermion. This annihilation process is not efficient enough to satisfy the new tight 
upper limit on the relic density of dark matter. One needs some additional contribution 
from the fi-coannihilation or Higgs exchange channels. In fact, the region at low Mq — M\i% 
is only allowed at large tan/3 because there one gets extra contribution from s-channel 
Higgs exchange as well as from coannihilation. As we have just discussed the low M —Mi/ 2 
region is also constrained from the Higgs mass limit at least for intermediate tan/3 and 
by the measurement of b — > 57 for tan (3 > 35. We also remark that the top quark mass 
plays an important role in determining the relic density constraint in mSUGRA models. In 
some cases the special tuning of parameters required to satisfy the relic density constraint 
can only be obtained for specific values of the top quark mass. For example two of the 
preferred regions, the focus point region and the heavy Higgs funnel regions, are very 
sensitive, through the renormalization group equations, to the value of the running top 
quarks mass. We first discuss in detail the case m t = 175 GeV before letting the mass 
vary within the la range. 



3.2.1 m t = 175 GeV 

For any value of tan/3, the regions consistent with the WMAP measurement correspond, 
for a fixed value of m t , to thin strips in parameter space. The WMAP allowed region 
corresponds to the very narrow area between the contours flh 2 = .094 and Qh 2 = .129, 
see for example Fig. El for m t = 175 GeV. We are then left mainly with three regions: 
fi-coannihilation, focus point and Higgs funnel. Although we will generally refer to the 
Higgs funnel region as the region where neutralino annihilation proceeds through the s- 
channel exchange of a heavy Higgs scalar, we note that there is also a small region where 
light neutralinos annihilate efficiently near the light Higgs resonance. For example, this 
can be seen as the narrow vertical strip along the LEP2 limit on the chargino mass near 
M 1/2 = 150 GeV, in Fig. 0for tan/3 = 10. 

In mSUGRA, the fi-coannihilation region designates the region at low Mq, where the 
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Figure 2: Allowed regions in the Mq — M1/2 plane for tan/5 = 5, 10, 35, 50, Aq = 0, fi > 
and mt = 175 GeV. Contours for = 111, 114GeV (green) Sa^ = 5, 10,30,60 x 10 -10 
(red) (top to bottom) and for flh 2 = .129(.094) (black/pink dots). The contour b — > 57 = 
2.25 x 10~ 4 (dash- dot/blue) is also displayed for tan (3 = 35,50. The region to the left of 
this contour is excluded. The LEP direct limits exclude the region delimitated by the blue 
vertical band. The regions between the Qh 2 = .094 contour and the theoretically excluded 
regions (yellow/light grey) have a value for the relic density below the WMAP bound. 
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lightest fi is the NLSP. As the relic density is very sensitive to the NLSP-LSP mass differ- 
ence, the proper mass degeneracy for a relic density consistent with WM AP corresponds 
to a very narrow strip in the M — M1/2 plane, see Fig. El Naturally this strip lies slightly 
above the forbidden region where < m^o. This strong dependence on the mass differ- 
ence also means that it is crucial to properly determine the pole masses from the running 
DR masses. In the coannihilation region, the main channels are Xi^i ~~ > r 7 an d T\f\ — > rf. 
As one moves towards degenerate f/xi, the relic density eventually drops below the 
WM AP range. Typically a value of the relic density within the WMAP range requires a 
mass difference AM^o ^5 — 15 GeV for m^o m (9(100) GeV while a smaller mass splitting 
is required for heavier neutralinos (near degeneracy is required for m^o = 400GeV). Fig. El 
shows the relic density versus the mass difference for tan (3 = 10, 35 as well as the mass 
difference as a function of the neutralino mass. At tan (3 = 35, one can afford a larger mass 
splitting as there is a larger contribution from the s-channel Higgs exchange to neutralino 
annihilation into fermions as well as an important contribution to x\^ ~ > The latter 
contains a term proportionnal to the sfermion mixing. The mass difference between the 
LSP and the sleptons of the first two generations can be much larger and ranges from 
10-20GeV for lower values of tan (5 up to 100 GeV when tan (3 > 35. There, the contri- 
bution from coannihilation processes with first generation sleptons become irrelevant 4 . 
Eventually, as the neutralino mass increases the coannihilation cross sections also become 
too small and one obtains an upper limit on the LSP mass. When combining this upper 
limit with the constraint imposed by the Higgs mass as well as from b — > sj, we find 
the range of allowed values for tan/5 = 10(35), 270(342)Ge\/ < M 1/2 < 1050(1250)GeV 
corresponding to neutralinos of 105(138)CreV < m^o < 448(540) GeV. 

The Higgs annihilation region includes both a light Higgs annihilation at low Mi/ 2 
and a new region which appears only for large tan (3. The former occurs only for the LSP 
mass very close to M^/2 since the light Higgs has a very narrow width. The latter appears 
at moderate Mq — Mu% values and corresponds to the annihilation of neutralinos via the 
s-channel exchange of a heavy Higgs. The main channels are X1X1 ~^ bb, rf corresponding 
to the preferred decay channels of the Heavy Higgs. The crucial physical parameter is 
the mass difference between Ma and 2m^o relative to the pseudoscalar width, hence the 
importance of including higher-order corrections for the width (7S|. Note that for very 
heavy scalars the width can be quite substantial, O(30)GeV. In FigEJ the mass difference 
AM^o = M4 — 2m^o relative to the width, T^, is displayed as function of M1/2 for 
tan f3 = 50. Here all constraints are applied. Models that fit within the WMAP range are 
roughly within AM^o < 2r^ except in the low M1/2 region. There the coannihilation 
region merges with the Higgs funnel region. Even though the coannihilation channels 
never dominate for these large values of tan/3 the contribution of the light sleptons to 
neutralino pair annihilation is sufficient to get an acceptable value for Qh 2 . In fact, nearly 
degenerate f/xi will lead to a very low value for the relic density. Due to the heavy Higgs 
annihilation channel, neutralinos as heavy as 650GeV can give reasonable values for the 
relic density. This type of model could be very difficult to hunt at the LHC as it features 
a rather heavy spectrum. Indeed at the tip of the Higgs funnel, one finds all squarks 
and the gluino in the 2.5 — 3.2 TeV range. Note that the position of the Higgs funnel is 
sensitive to the standard model input value used in the RGE codes, in particular mb(mb) 
[HS] as well as the top quark mass. Indeed the soft scalar masses in the Higgs potential, 

Note that the small mass difference between the stau and the LSP raises the question of detectability 
of a nearly degenerate f lepton as well as the feasability of a precise measurement of its mass at a collider. 
For recent studies within the context of a linear collider see PHI EH • 
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Figure 3: a) flh 2 as a function of the mass difference AM^o in the coannihilation region 
for tan/3 = 10(black dots) and tan/3 = 35(green crosses), b) For models where Qh 2 is 
within the WMAP range, m^o as a function of the mass difference AM^o for tan (5 — 10 
(black dots) and tan/3 = 35 (green crosses) and as a function of the mass difference 
AM^o for tan (3 = 10 (blue dots). 



(m^ ) , from which one evaluates the mass of the pseudoscalar at the EWSB scale are 
very sensitive to the bottom (top) Yukawa's [To] . 

The last allowed region is the focus point region, this region is found at high values 
of M where the value of /1 drops rapidly. When /1 < M 1; M 2 , the LSP has a significant 
Higgsino fraction. Furthermore the NLSP (x\ and/or xf) are also mostly Higgsino and are 
not much heavier than the LSP. Thus the neutralinos/charginos coannihilation channels 
are favoured. The position of the focus point region is very sensitive to the value of the top 
quark mass that enters the RGE and also differs for different RGE codes [TBI IKJ1 E2] • For 
the case of SOFTSUSY used here, rather low values of the top quark mass are necessary 
to reach the region where /1 drops rapidly, at least for intermediate values of tan /3. As 
displayed in Fig. [21 for m t = 175 GeV, one finds a cosmologically allowed region at large 
M only for large tan/3. In the focus point region, the main annihilation channels are 
X1X1 ~^ W + W~ , ZZ. Coannihilation channels such as X1X1 ~~ > QQ* can a l so help reduce 
the relic density. One also finds an important contribution from X1X1 ~~ * tt as soon as this 
channel becomes kinematically accessible. This occurs, in the focus point region, when 
Mo ~ 4.5 — 5TeV, Fig. EJl. Even though the stops are heavy (777^ w 2.5TeV), a sufficiently 
large annihilation cross section results from the Z exchange contribution. For the gauge 
bosons annihilation channels, values consistent with WMAP for Qh 2 require not so large 
a value for /x This implies that the chargino and lightest neutralinos should be rather 
light. For example, the chargino and neutralino are confined to have 777-+ < 375GeV 
for tan/3 = 50, 777 t = 175 GeV. The chargino and neutralino can be nearly degenerate 
O(10)GeV but this occurs mainly when the relic density falls below the WMAP range. For 
flh 2 ~ .1 typical mass differences are rather around 50GeV. Note that models in the focus 
point region feature a very heavy sfermion sector, the LHC will have little opportunity 
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to discover the squarks, only the gluino could be accessible. At the same time the rather 
light Higgsino sector could be probed at a linear collider, extending the reach of the LHC 
[HI] • Furthermore these models are also interesting for future direct dark matter detection 
experiments as will be discussed in Section 



3.2.2 m t + 175 GeV 

The position of the allowed region at large Mo is the main qualitative difference between 
mSUGRA models with a lower input values for the top quark mass while the more strin- 
gent constraint originating from the light Higgs mass shifts the range of allowed values 
for sparticle masses in the low Mo region. For example, when m t = 170GeV, one finds 
a focus point region even at tan ft = 10, Fig. At large values of tan ft, the focus point 
region becomes much wider. Indeed all points between the contour delimiting the region 
forbidden by the requirement of electroweak symmetry breaking (EWSB) and Qh 2 = .094 
have a value for the relic density below the WMAP range, indicating eventually some 
other dark matter component. The shape of the allowed region at tan/5 = 50 differs 
significantly from the case with heavier top quark mass. First, one tends to get a lighter 
heavy Higgs, extending the Higgs funnel region towards higher values of Mo, at the same 
time shifts in \i are found pushing the focus point towards lower values of Mo. This is 
solely an effect of the RGE dependence of the pseudoscalar mass and of \i on the top 
quark mass. As concerns the relic density constraint the requirement remains the same, 
only models where the neutralino lies near Ma/2, roughly AM/T^ ~ 2 are allowed. The 
Higgs funnel and focus point regions merge for Mu2 ~ 1 TeV. Furthermore, because of 
the additional contribution from Higgs exchange, one finds acceptable values for the relic 
density in the low Mo — M1/2 region. In the end, a large fraction of parameter space is 
allowed by the relic density constraint. Of course direct measurements such as b — > 57 
play an important role here particularly in the low Mq — M]/ 2 region. 
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Figure 5: Same as in Fig. 2 for m t = 170 GeF 
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Conversely a heavier top quark, say m t = 179GeV, means a shift of the focus point 
region towards higher values of M [S3]. We have stopped our scans at M = 5 TeV, 
as we found it rather difficult to obtain a solution to the RGE that converged rapidly 
for very large Mq 5 . The relaxed bound from the Higgs mass also favours the region of 
parameter space where not so heavy neutralinos/charginos and sleptons can be found, here 
m^o > 90 GeV for tan/3 = 10. Finally the Higgs funnel moves closer to the coannihilation 
region, indeed the heavy Higgs masses are shifted upwards with the increase of the top 
Yukawa making it increasingly difficult to have m^o pa Ma/2. Then for tan/3 = 50 the 
allowed region only extends slightly above the stau coannihilation region. 

mSUGRA, the case A ^ 

For simplicity we have up to now fixed the value of the trilinear coupling to Aq = 0. 
Although the constraints are not extremely sensitive to this parameter yet as one varies 
Aq over the full range \Aq\ < 2TeV, one sees an impact on the allowed regions in the 
Mq — Mi/2 plane. For one, a non-zero value for the trilinear coupling shifts the physically 
allowed region where the neutralino is the LSP. First, a region with tachyons appears at 
low Mi/ 2 . Furthermore, the stau mass is lowered since the trilinear mixing affects the 
running of the t^r masses as well as the running of fi which determines the mixing in the 
stau sector. Indeed the mixing is proportionnal to A T + ft tan (3 so is usually dominated by 
the term \x tan (3. The boundary of the coannihilation region then moves towards higher 
values of Mq for non-zero values of Aq. This is particularly obvious at large tan /3, Fig.0|8j 

First consider models with A = 2000 GeV and m t = 175 GeV. One finds that 
consistency with WMAP is limited to a region just at the boundary of the theoretically 
excluded region for tan/? = 10, see Fig. One finds a new feature in these scenarios, a 
region where the t\ is the NLSP. This region which occurs at moderate values of Mq lies 
just to the right of the region where one finds tachyonic solutions, i.e. for 200 < Mq < 
650 GeV. The main annihilation for neutralinos is into fermion pairs and coannihilation 
channels such as Xi^i — ¥ tg, W + b contribute also significantly. Note that the contribution 
from Xi^i ~> W + b can give the dominant contribution to the thermally averaged cross 
section exceeding the QCD channels with gluon emission 6 . Then one can satisfy the relic 
density constraint because of the coannihilation processes involving stop quarks. However 
this stop coannihilation region is not allowed by the light Higgs mass constraint. This 
constraint is more severe than in models where the trilinear coupling vanishes, only the 
stau coannihilation region remains once the limit is imposed. The allowed values 
for neutralinos masses in the coannihilation region lie in the range 170(126) < m^o < 
440(700)GeV for tan/3 = 10(50). Of course one can have models with an even heavier 
LSP since for tan/3 = 50, both the focus point region as well as the Higgs funnel are 
present. These regions are shifted as compared to the A = case, in particular the heavy 
Higgs masses are shifted down making for a more efficient annihilation of neutralinos via 
a heavy Higgs exchange. Then neutralinos as heavy as m^o pa 740 GeV are allowed. 

For negative values of A , say A = — ITeV, the constraints from the Higgs mass are 
relaxed while the b — > s'j is more important as was discussed in the previous section. The 
relic density is the main constraint on models with low tan (3. Again one finds a region 
where the t\ is the LSP at low Mq — Mu% values. Just when the t\ becomes the NLSP, 

5 The predictions of the different codes for the SUSY spectrum calculation can be very different in this 
region PJ- For example, with 1SAJET the authors of [EE] find the focus point region all the way up to 
M = lOTeV. 

6 This contribution is not included in |8fi| . 
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values of the relic density consistent with WMAP are found, see FiglHlfor tan/5 = 10. 
However this region is barely compatible with the Higgs mass limit and does not pass the 
b — > S7 constraint. One also finds as usual a thin f\ coannihilation region. Furthermore, 
with the large negative mixing which tends to increase the heavy Higgs masses, the Higgs 
funnel is reduced to a tiny region even at tan/3 = 50, Fig.|H| Note also that the constraint 
from B s — > starts to restrict the mSUGRA parameter space although it is not yet 

competitive with other direct constraints. For larger values of the mixing parameter, e.g. 
Aq = —1700 GeV, it is possible to find models that satisfy all constraints and where the 
tx is the NLSP jEZj. 



3.3 x 2 fit 

So far we have discussed the effect of each type of constraint taken individually as well 
as their sensitivity on the top quark mass. To get a better estimate of the impact of 
the uncertainty on the input parameters, we have also performed a y 2 fit to the four 
observables: Qh 2 ,b —> 57, rrih and mt for a fixed value of tan/3 and Aq. We use m t = 
174.3 ±5.1 GeV. We show in Fig. ©the allowed area at 68% and 95% C.L. One finds that 
the favoured regions are the coannihilation as well as the focus point region for tan j3 = 10. 
The area of parameter space allowed at 95%C.L. after allowing for the m t uncertainty now 
seems very wide especially at M > 2 TeV. This reflects the high sensitivity of the RGE to 
the value of the top quark mass. For tan (3 = 50, the fi-coannihilation, Higgs funnel and 
focus point region are not easily distinguishable and a large fraction of parameter space 
is allowed. Note that both the Higgs funnel and focus point regions are very sensitive to 
the value of the top quark mass as we have seen in Fig. 
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and B s — > fi + fi~ = 2 x 10 -7 (dots) for tan(3 = 50. The exclusion regions are labelled in 
the same manner as in Fig. [B 
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In summary, in mSUGRA models a combination of the direct collider limits and relic 
density constraints, restricts the allowed models in very specific regions of parameter 
space. The exact location of these regions in the M — Mu2 plane is however sensitive 
to theoretical uncertainties in the evaluation of the spectrum. Nevertheless the predic- 
tions at the EW scale are rather contrived: either nearly degenerate f /neutralinos, light 
higgsino-type neutralinos or at large tan/3, neutralinos that are near Ma/2. The results 
we have presented agree qualitatively with other published results, however differences in 
the position of the heavy Higgs funnel as well as the focus point region are found jH IHH] • 
To a large extent these differences can be traced back to differences in the evaluation of 
the supersymmetric spectrum at the EW scale. These issues were brought up in Ref. [3T] 
and more details will appear in [32] . 

4 Non-universal gauginos 

In models where one relaxes the gaugino universality conditions the allowed regions in 
parameter space can change drastically |89 ] 190 ] l9"T] . For one, the wino content of the LSP 
can be significantly increased, favouring the annihilation cross section into pairs of gauge 
bosons. This occurs for example for M[/M 2 > 1 at the GUT scale. Then reasonable 
values for the relic density can be obtained even for heavy sleptons. The annihilation 
cross section into pairs of W tends to be more efficient than in lepton pairs [22] in fact 
so efficient that the relic density often falls below the WMAP range. Furthermore, the 
coannihilation with heavier charginos and neutralinos, also the most efficient ones, are 
much more frequent in this type of model. Second, as in the universal case, one finds 
significant areas of parameter space where neutralino annihilation proceeds through a s- 
channel Higgs exchange. However these regions do not necessarily only appear at large 
tan/3. Allowing M^/M 2 < 1 at the GUT scale reduces the value of the heavy Higgs 
masses without much effect on the neutralino mass. This type of models also predicts 
smaller values for fi which means a LSP with a significant Higgsino component and more 
efficient annihilation channels without requiring high tan/3 Alternatively, reducing 
M[/M 2 allows annihilation of neutralino through a light Higgs exchange. Finally non- 
universality in the gaugino sector can also affect the predictions for the sfermion masses. 
In particular reducing the parameter M' z for given M[, M' 2 values, produces much lighter 
squarks. Coannihilation channels with squarks then are more likely especially when there 
is a large mixing in the squark sector. Altogether one finds the same mechanisms as in 
mSUGRA for getting sufficient annihilation/coannihilation of neutralinos but the allowed 
region in the parameter space of the MSSM shifts significantly and no longer requires a fine 
tuning of parameters. The direct constraints are of course also modified in non-universal 
models. This will be detailed next and illustrated with a few typical case studies. Here 
we only consider the case m t = 175GeV. 

4.1 Mg > M' 2 , the case M3 = 2M' 2 = 2M[ 

By changing only M3 while keeping M 2 = M[ at the GUT scale means that at the weak 
scale one retains the usual relation Mi ~ 0.4M 2 . Futhermore the parameter \x tends to 
be larger than in mSUGRA, see Eq. Therefore the Xi is nearly a pure bino. The relic 
density constraint then becomes more severe than in the universal case since annihilation 
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of bino LSP's into light fermions, although dominant, is not sufficient to bring Qh 2 in the 
WMAP range. The pseudoscalar mass, which is driven by M 3 , is also much heavier than 
in the universal case so annihilation through the heavy Higgs cannot take place. The 
focus point region is pushed towards higher values of Mq even at tan/? = 50. Then the 
only allowed region is the fi/xl coannihilation region. This region is always present since 
the relation between Mq and Mi, M 2 which to a large extent determines m^o, has not 
changed. The allowed narrow band extends beyond M1/2 = 2TeV, see so that a 

rather heavy SUSY spectrum is possible. In the coannihilation region one finds roughly 
AM f ^o ^5 — 10 GeV for tan/5 = 10, comparable to the values found for the universal 
case. For larger values of tan (3, the increased contribution of many channels, in particular 
those with a light Higgs, x^f — > rh, ff — > hh means that larger mass differences can be 
compatible with WMAP. 

On the other hand the constraints from colliders are not as severe. As the Higgs 
mass is driven by the squark mass which increases with M3, lower values of M\i 2 = M' 2 
should be allowed. The upper bound on the Higgs mass is indeed easily satisfied for 
tan f3 = 10 as can be seen in FiglTOl Allowed regions are found even at tan f3 = 5 . 
There for rrih > 111 GeV one recovers a small allowed region with M\/ 2 > 250 GeV 
and light sleptons. The constraints from precision measurements are easily satisfied in 
this model since heavy squarks also imply that b — > 57 is close to the standard model 
prediction. Only a tiny region is excluded for tan/3 = 50. In this type of model, for 
example when My 2 > 800 GeV, the heavy squarks make it difficult for the LHC to 
discover supersymmetry. On the other hand one recovers part of the low Mu 2 region that 
is favourable for collider searches of charginos and sfermions. 




M 1/2 (GeV) M 1/2 (GeV) 

Figure 10: Same as in Fig. \^for fi > 0, Aq = 0, m t = 175 and M3 = 2M' 2 . Only contours 
for nih = 114 GeV (green) are relevant. 



4.2 M3 < M' 2 , the case M3 = M' 2 /2 

Decreasing the value of M3 while keeping M 2 = M[ at the GUT scale means that the 
relic density constraint can be satisfied more easily. First, one finds a neutralino-LSP 
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with higher Higgsino component, although the LSP remains mainly a bino. Indeed /i is 
typically smaller than in unified models, EqE3 More importantly one feels the presence 
of the heavy Higgs pole even for moderate values of tan/?. As we have discussed in 
Section ^ solutions for Ma ~ 2m^o can be found even for tan/5 = 10, Eq. |H1 Neutralino 
annihilation into fermions is then dominated by the s-channel scalar exchange and the 
main channels are X1X1 — > ti, bb, rf, the former being dominant at lower values of tan /3 
only. Note that for tan/3 = 10 the coupling of the pseudoscalar to b quarks is already 
2 to 3 times larger than to t quarks the ratio ^| ~ ^ tan (3 2 and this mode starts to 
dominate. The region that in the universal case was called the f coannihilation region 
now becomes a Higgs/coannihilation one and includes the low M — M\/2 region, similar 
to what was observed previously for large values of tan/3. In fact as one moves closer 
to the f LSP region and coannihilation channels start to contribute significantly to the 
effective cross section, the relic density is already much below the WMAP range. This 
corresponds in Fig^^to the region between the f LSP line and the contour Qh 2 = .094. 
This means that values of Qh 2 within the range of WMAP are no longer associated 
with nearly degenerate f — Xi- As in mSUGRA, as one moves towards larger values of 
Mi/2 the stau coannihilation cross sections become much smaller and even models with 
completely degenerate stau/neutralino cannot meet the upper limit of WMAP on the relic 
density. However a new Higgs funnel region appears, it corresponds again to s-channel 
Higgs resonance but this time for the coannihilation processes X1X2J X1X3 ~^ tt, bb, rf as 
well as charged Higgs exchange in X1X1 ~ * tb. The onset of the charged Higgs exchange 
contribution is easily visible as a kink in the contours of constant relic density in Fig. 
ITTk-b. The focus point region is also more important even for tan/3 = 10. 

As one increases tan/3, the heavy Higgs decay into b-quarks completely dominates. 
The allowed region in the Mo — M1/2 plane correspond mainly to annihilation into bb, rf. 
One does not find a region where f-coannihilation dominates, although f — Xi can be 
quite close in mass, this occurs only when one is near the Higgs resonance. Then the relic 
density is much below the WMAP range. For values of M1/2 > HOOGeV, only annihilation 
or coannihilation with heavier neutralinos/charginos into fermion pairs through a Higgs 
exchange take place. Since the value of \i is smaller than in the universal case, there is also 
a more important contribution from the annihilation into gauge bosons channel and the 
focus point region is much wider for a fixed top quark mass. Eventually the focus point 
and Higgs annihilation region merge, for example for tan/3 = 35 and for M1/2 ~ 2TeV 
the entire range of Md is below the WMAP upper bound. For tan/3 = 50, the merging 
of the two regions is complete and a large fraction of the parameter space is allowed. As 
for other values of tan/3, annihilation into bb pairs dominates for low values of Mo and 
annihilation (or coannihilation) into gauge boson pairs eventually takes over as one moves 
close to the region where EWSB cannot take place. We remark that even with a moderate 
trilinear mixing, say A = — 600GeV, one finds regions where the t\ is the NLSP. 

As one would expect, when M3 < M' 2 , the Higgs mass constraint becomes more severe 
as the squark masses are lower than in the universal case for a fixed value of Mu% = 
M' 2 . For example, for tan/3 = 10, only My 2 > 460Ge\^ is allowed in the coannihilation 
region. The lighter squarks relative to the universal case imply a larger deviation from 
the standard model prediction for b — > sj as well. The exclusion region in the left-hand 
corner of the M — M1/2 plane where the b — > 57 branching ratio drops below the allowed 
range becomes more important. As usual this constraint takes over the one coming from 
the Higgs mass as one increases tan/3. For Sa^, the main difference from the universal 
case is that the value of /1 are modified by the RGE's. Nevertheless the allowed region 
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is very similar to the one obtained in the universal case, the contours are only slightly 
shifted towards higher M\i 2 values. Note that here again requiring a strict deviation from 
the standard model, would rule out most of the region at large M — M\i 2 . Combining 
this with other direct constraints would leave only a small region of parameter space for 
tan (3 < 10 as the direct constraints are more effective in this class of model. 

In these models typically squarks are lighter than in mSUGRA, one then expects 
a good discovery potential for the LHC However as one easily finds a region where 
m^o pa Ma/2, even for very heavy neutralinos one can find reasonable values for the relic 
density. It is possible then to obtain a very heavy supersymmetric spectrum. For example 
for Mq = 4TeV, M\/ 2 = 3TeV, the spectrum consists of heavy squarks, rriq = 3 — 4TeV, 
neutralinos and charginos, m^o pa m-+ ~ 820GeV and also heavy Higgses m# pa itih+ ~ 
rriA ~ 1.5TeV, a difficult task for discovering supersymmetry at colliders. The only new 
particle that could be reached at a collider would be the light Higgs. In these models even 
Safj, does not differ significantly from zero. Only an unambiguous evidence for a non-zero 
contribution to 5a^ would restrict this class of models with heavy sparticles. 

4.3 The general case M' Z ^M' 2 = M[ 

There is a smooth transition from the region with universality where Qh 2 is in general too 
large to the region M' z < M 2 where there is not enough relic density. This is illustrated 
in Fig. ^1 showing the value of Qh 2 as function of r 32 for three typical values of M = 
300, 1500, 3000 GeV. For these values of M one crosses the coannihilation, Higgs exchange 
and focus point region respectively for m t = 175 GeV. Clearly it is much easier to satisfy 
the relic density upper bound in models where Mg < M 2 . As just discussed this is the 
effect of the heavy Higgs pole and of the more important Higgsino component of the 
bino LSP. The Higgs funnel can be found even for Mq = 300 GeV as long as r 32 < -5. 
Furthermore for tan /3 = 50, and Mq = 3000 GeV, nearly all models are below the WMAP 
range when r 32 < -4, this is a combined effect of a smaller value for /i and Ma as compared 
to the universal case. 

In conclusion, relaxing the unification condition on M 3 /M 2 while keeping M[ = M 2 
helps recover some values for low tan/? because of either a weaker Higgs mass limit or 
possible annihilation through a heavy Higgs resonance. This was illustrated for tan (3 = 5. 
Furthermore, the case M 3 < M 2 satisfies more easily the relic density constraints for any 
values of tan (3. In particular for large tan f3 where the relic density constraint is satisfied 
in a large portion of the available parameter space. A significant fraction of the r 32 < 1 
models then predict squarks within the range accessible by the LHC. On the other hand 
because of annihilation/coannihilation through the heavy Higgses, scenarios with heavy 
neutralinos as well as heavy sfermions can satisfy all constraints while predicting few 
signals at colliders. In these models a value of da^ compatible with the standard model 
is predicted, evidence of a deviation might be the best way to rule them out. 

4.4 The case M' 2 = M3 

We now restore the unification condition M 3 = M 2 while allowing M[ to differ from 
M 2 at the GUT scale. This will directly affect the wino/bino nature of the LSP. If the 
LSP has a higher wino component, which occurs in models where M\ > M 2 at the weak 
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Figure 11: Same as Figure 2 for \i > 0, A = 0, m t = 175 GeV and Mg = M' 2 
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Figure 12: Vth 2 vs r 32 for tan/3 = 10, M = 300, 150QGeV and tan/3 = 50, M = 
1500, 3000 GeV. In all cases, fi > 0, A = 0, mt = 175 GeV and the constraints from 
b — > sj, nih as well as LEP direct limits are imposed. 



scale, annihilation into W pairs becomes important especially since LEP constraints forces 
m-jj.o > M\y. Furthermore this process is much more efficient than the annihilation into 
fermion pairs that dominated for a bino LSP. A wino LSP also opens up the coannihilation 
channels with other gauginos. Then one expects the relic density constraint to be easily 
satisfied in models where M[ > M' 2 . For fixed values of M , we scanned over models 
with .2 < ri2 < 5 and tan/3 = 10,50. Figll 31 clearly shows the dramatic effect on Qh 2 of 
increasing ri 2 , varying from almost always too large relic density when ri 2 < 1 to nearly 
always too small when ri 2 » 1 passing by an intermediate region where flh 2 lies within 
the desired range. This value preferred by WMAP is around r\i = 1.8 for tan/3 = 10 
and corresponds to Mi ~ M 2 at the electroweak scale. Thus the lightest neutralino is 
a mixed bino/wino. We will study this case in more details below. Note that we have 
stopped our scans at M 2 = 2TeV so that we have not allowed very heavy neutralinos. 
Had we extended the parameter space it would have been possible to find models with 
very heavy neutralinos that fell within the range of WMAP even for r 12 > 2. For higher 
values of tan/3 roughly the same behaviour is observed, only then a wider selection of 
models passes the constraint. For example, for tan/3 = 50 and Mq = 1500GeV, many 
models satisfy the relic density constraint, this reflects the presence of the Higgs funnel 
which appears already for ri 2 « 1. 
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Figure 13: Relic density vs ryi in models with fi > 0, Aq = 0, mt = 175 GeV for a-c) 
tan/3 = 10 and M = 300, 1500, 3000 GeV d) tan/3 = 50, M = 1500 GeV. The horizontal 
lines denote the 2a region of WMAP. 



4.5 M[ > M' 2 , the case M[ = = 1.8M£ 

We first consider the special case M[ = I.8M2. This type of model can arise from 
string theory. For example in models where gaugino masses are loop suppressed and non- 
universal one loop contribution become dominant, models have been constructed where 
at the one- loop level the mass M[ exceeds M' 2 and M' z ESI- We recall that at the 
electroweak scale this model corresponds to M\ m M 2 . The lightest neutralino while 
still dominantly gaugino has a large wino component and is typically heavier than in 
mSUGRA. Thus the increased importance of the neutralino annnihilation into W + W~ . 
Furthermore the LSP can be nearly degenerate with the chargino and xl- This favours the 
gaugino coannihilation channels. The main channels are into quarks and gauge bosons, 
XiXii X1X2? X1X1 ~* W + W~ and XiXi ~^ qq,ZW. The allowed region is not anymore 
restricted to the low M region as even when squarks and sleptons are very heavy so that 
the fermion annihilation channels becomes small the annihilation into gauge bosons still 
occurs and has a large enough rate. In fact this mode is so efficient that this type of model 
often leads to a value for the relic density that is too low. For example for tan/3 = 10, in 
Fig. El the whole region to the left of the flh 2 = .094 contour has a relic density below 
the WMAP range. Agreement with WMAP then implies a rather heavy LSP, for example 
m,o w 600 GeV for tan 3 = 10. 

For neutralinos up to roughly M1/2 = 700 GeV, the coannihilation with the f\ implies 
a relic density below the WMAP range. Above this value one does not find an important 
coannihilation region. The region between the stau LSP border and the Qh 2 = .129 
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contour features too high a value for the relic density. There the NLSP is a slepton but 
slepton coannihilation exchange is not efficient enough for very heavy LSP as we have seen 
in the mSUGRA case. Since, for a given M1/2 the LSP is much heavier than in mSUGRA, 
the scale is sligthly misleading. As one increases Mq and the mass of the chargino becomes 
more degenerate with the mass of the neutralino then the coannihilation with chargino 
becomes significant and there is a region that satisfies WMAP, in this region one is close 
to the heavy Higgs pole. The kink at large Mi/ 2 occurs when the coannihilation process 
XiXf ~^ tb lies near the charged Higgs resonance. For large tan (3, the heavy Higgs pole 
moves towards lower M\m and merges with the bulk region as well as the focus point 
region. 

The constraints on this type of models from the Higgs mass as well as from precision 
measurements are not significantly different than in the universal model. For example, 
the allowed region starts roughly from M1/2 > 300(450)GeV for tan/3 = 10(50). The first 
bound is set from the Higgs mass, the second from b — > sj. Note that the Higgs mass 
and the b — > 37 constrain only models where the relic density is below WMAP. 

For collider searches, the spectrum of sparticles as concerns the coloured sector is not 
dramatically different then in mSUGRA. However since the WMAP allowed regions are 
completely different and in particular include a Higgs funnel at small tan /3, models with 
very heavy squarks are perfectly acceptable. For example, for M = M\/i = 2 TeV and 
tan/3 = 10, all squarks are in the 3-4 TeV range. These are beyond the reach of LHC. 
Furthermore, the LSP tends to be much heavier than in mSUGRA, more importantly, 
there is a small mass difference between the chargino and the neutralino which can effect 
searches. 

Had we increased further the ratio 7*12 we would have found a LSP with an even 
larger wino component. This would increase the values for the LSP masses compatible 
with WMAP. In this way one can easily reconcile the relic density constraint even with 
a very heavy SUSY spectrum. For example if r 12 = 2.5 the LSP as well as the whole 
super symmetric spectrum lies above 1.6 TeV. 

4.6 M[<M' 2 = M3 

In that case, the neutralino can be very light as the constraint on the chargino implies a 
lighter neutralino than in mSUGRA. A general analysis of MSSM models with light neu- 
tralinos and including all constraints has shown that neutralinos as light as 6 GeV could 
be allowed provided the heavy Higgses were rather light 12Z] • In SUGRA type models 
with non-universal gaugino masses however the heavy Higgs mass is not an arbitrary pa- 
rameter but is deduced from the RGE. Typically one finds a rather heavy scalar, so one 
cannot appeal to s-channel annihilation through a heavy Higgs resonance to bring the relic 
density in the desired range, mainly the light Higgs exchange is available. As an example 
take M' x jM' 2 = 1/3, the main annihilation channel is into a pair of fermions. One finds, in 
addition to the f coannihilation region, a region allowed by the relic density constraint, 
that corresponds to values of Mi/ 2 where the neutralino can annihilate through a light 
Higgs pole. This region is compatible with the light Higgs mass limit which is more or less 
the same as in universal models. For large tan j3, as in universal models, annihilation into 
fermions receives a contribution from the heavy Higgs exchange but this is not sufficient 
to allow very light neutralinos. The relic density constrains neutralinos to be heavier than 
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m^o & 30GeV. This type of model offers good prospect for direct detection of the LSP 
as well as for supersymmetry discovery both at the LHC and at the linear collider as was 
discussed in j2HJ 123 ■ 

4.7 Other non- universal cases 

Finally there is another class of models that feature non-universal gaugino masses, models 
with M[ = M3 7^ M 2 EU inn] ■ We do not discuss these models here as the main impact 
of relaxing universality has already been illustrated in the cases studied above in details. 
First we expect that when M[ = M3 > M 2 the relic density would be easily satisfied 
because of the large wino component. Of course, as compared to the case where M3 = M 2 
discussed above, shifts in the value of rriA and /x will result from an increase in M' z . Second 
we expect models with M[ = M3 < M 2 , to be rather similar to the M3 < M 2 models 
since the crucial parameters for the evaluation of the relic density, such as \x and M4 are 
determined by M' z via the renormalization group equations. 

5 Prospects for direct detection 

We have just shown that cosmological measurements provide strong constraints on many 
super symmetric models. Here we discuss the potential of direct detection searches for the 
LSP both in the mSUGRA model and non-universal models. 

Several experiments have placed limits on the cross section for the spin independent 
neutralino proton scattering, typically at the level of 10 _6 pb [33]. There is also one report 
of a signal by DAMA|35]. In the near future results form several upgraded detectors are 
expected. In fact, the first results from CDMS2 have been published recently jSE] showing 
already significant improvement on the upper limit. Finally ton-scale detectors are also 
planned, such as ZeplinIV[37j, Genius [38 j and Xenon[3^. These large scale detector 
will reach a sensitivity between 10~ 9 — 10~ 10 pb for neutralinos around lOOGeV. Typically 
detectors have a reduced sensitivity as one increases the neutralino mass. In parallel, 
detectors that are also sensitive to the spin dependent part of the neutralino proton cross 
section have obtained their first results, NaiadpEO]. SimplejH] and PICASSO 42 . At the 
moment these detectors are much less precise and typically limits at the level of a few 
pb's are given. Upgrades of these detectors are planned and eventually a ton-size detector 
such as PICASSO could reach a sensitivity between — 10~ 4 pb. 

We follow the usual practice and use the cross section for annihilation of neutralino on 
proton, cr^o p , to compare the direct detection potential in a given model. We do not rescale 
the cross section to take into account the possibility that the local density of neutralino 
could be lower than expected. However one has to keep this point in mind especially 
for the models where the relic density is very low. The spin-independent cross section 
is dominated by the Higgs exchange and squark exchange diagrams whereas the spin 
dependent one receives contribution from the Z exchange and squark exchange diagrams. 
One therefore expects to find the largest cross sections, for both spin-dependent and spin- 
independent processes, in the low M region where light squarks are found as well as in 
the focus point region of mSUGRA models |07J where the coupling of neutralinos to both 
the light Higgs and the Z are enhanced because of the mixed bino/Higgsino nature of the 
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neutralino. Note that in the focus point region, the large coupling to the Z/Higgs implies 
a large annihilation cross sections of neutralinos as well as large neutralino proton cross 
sections. 

5.1 mSUGRA 

We have performed scans over the parameter space of the mSUGRA model with \x > 0, 
A = 0, M < 5TeV and My 2 < 2TeV for tan/3 = 10 and tan/3 = 50 as in previous 
sections. We first consider a top quark mass m t = 175GeV. While the prediction for 
the neutralino proton cross sections vary over several orders of magnitude over the full 
parameter space, after imposing the WMAP upper bound one distinguishes easily the stau 
coannihilation region, the light Higgs annihilation region and the focus point region. In 
the latter case the cross sections are typically larger despite the suppressed contribution 
of the squark exchange. The reason is the enhanced coupling of the neutralino to the 
Higgs due to the mixed Higgsino/gaugino nature of the neutralino. Note that the same 
mechanism that made for efficient neutralino annihilation also induces large values for the 
spin independent cross section. As discussed previously, the focus point region is found 
only at large tan/3 in mSUGRA models for m t = 175 GeV. In Fig. fTHI the focus point 
region for tan/3 = 50 corresponds to the set of points clustered around m^o = 300 GeV 
for which a,si ~ 1 - 3 x 10~ 8 pb. This region can then be partly probed by EdelweissII 
and will be completely covered by future detectors such as ZeplinlV. In general, one also 
gets the contribution of the squark exchange diagram. Thus large cross sections are also 
expected in the low M Q region of parameter space where one finds the lighter squarks, this 
basically means the stau coannihilation region. Furthermore, the spin-independent cross 
section should increase for lighter neutralinos as they are found in the low Mq—Mi/2 region 
where in addition to having light squarks the coupling of the neutralino to the Higgs is 
large as well. Note however that this is precisely where the collider constraints come into 
play. After imposing all constraints we find the maximum value a$i < 3 x 10 _9 (10 _8 )pb 
for tan/3 = 10(50). These can only be partially probed by future large scale detectors. 
For example, Xenon will probe models for which m^o < 250(360) GeV. 

At large tan/3 there is also a region where neutralinos annihilate via a heavy Higgs 
exchange, in this region, neutralinos are heavy as well as the squarks and the heavy 
Higgs which lie in the TeV range. One therefore does not expect a large cross section 
despite the enhanced coupling of the Higgs to bottom quarks. Indeed this region features 
cross sections (9(10 _10 )pb and shows up as the continuation of the coannihilation region 
for m^o ~ 500 GeV in Fig. ED Finally the region at intermediate M where it is still 
possible to have annihilation of neutralinos via the exchange of a light Higgs scalar predict 
values within the range 10~ 10 (10~ 9 ) < a SI < 10~ 9 (Kr 8 ) pb for tan/3 = 10(50). The 
enhancement at large tan (3 is due to the heavy Higgs contribution. While the heavy Higgs 
exchange is completely negligible at intermediate values of tan /3 it eventually dominates 
the light Higgs despite the mass suppression factor. This is due partly to a decrease of 
the mass of the Heavy Higgs in models with large tan/3 but mainly one benefits from 
the enhanced couplings of the heavy Higgs to bottom quarks. This region is accessible to 
future large scale detectors such as Xenon j^Sj- 

The spin dependent cross section has roughly the same behaviour as the spin inde- 
pendent one. Indeed the neutralino LSP couplings to the Z is also enhanced in the focus 
point region because of the larger Higgsino component. This region corresponds to the 
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sest of points clustered around cjsd ~ 10~ 4 pb in FiglToTi and can be probed entirely by 
future large scale detectors such as PICASS0[42j. In the coannihilation region and the 
heavy Higgs annihilation region, the cross section decreases with an increase of the neu- 
tralino mass. As above this is related to the additional suppression of the heavy squark 
exchange. The predictions for the spin dependent cross sections still lie orders of mag- 
nitude below present and future limits especially in the heavy Higgs annihilation region 
at tan/? = 50 where <jsd ~ 10~ 8 pb. Finally the allowed models with light neutralinos 
for which a good relic density is found since neutralino annihilates via the light Higgs 
have also an enhanced coupling to the Z. Here one finds rather favourable cross sections, 
a si ~ — 1CT 5 . These models can be probed by a ton-scale detector. 

The neutralino proton cross sections should not depend very much on the value of 
the top quark mass for a given set of MSSM parameters. However the allowed parameter 
space and the predictions for the MSSM parameters within the mSUGRA model differs 
significantly as discussed in sectional If we increase the top quark mass, m t = 179 GeV, 
the focus point region will be shifted to higher values of M , out of the range scanned here. 
Furthermore the heavy Higgs annihilation region almost disappeared. We then recover 
the results for the spin independent cross section for mt = 175 GeV in the coannihilation 
and Higgs annihilation regions. However since the Higgs mass constraint is relaxed one 
can reach higher cross sections for light neutralinos, for example asi ~ 5 x 10 _8 pb for 
tan/? = 10. For tan/3 = 50 the maximum cross section in the coannihilation region is 
similar to the case of the lighter top quark since the constraints from colliders are much 
less dependent on the top quark mass. 

We also discuss the case of a rather light top quark since WMAP allows a more 
important region of parameter space and we can better illustrate the generic behaviour 
of the neutralino proton cross section in the different regions of mSUGRA models. The 
results of a scan over models with m t = 170GeV are displayed in Fig.E3 The focus point 
region with a Higgsino LSP is allowed even for tan (3 = 10, and as above that is where one 
finds the largest cross sections, asi ~ 10~ 8 — 10~ 7 pb over the allowed range of neutralino 
masses. Much lower cross sections are found in the coannihilation region. For tan f3 = 50, 
the allowed region covers most of the parameter space (see FigEb)- Nevertheless one 
can roughly distinguish the focus point, coannihilation and Higgs annihilation regions in 
Figllfib. The largest cross sections are again found in the focus point region where for 
light neutralinos asi can reach 10 _7 pb. The maximum value for the spin-independent 
cross section increases with tan f3 since one finds a neutralino LSP with a higher Higgsino 
content. In the coannihilation region which shows up as the lower branch in Fig. 116b. the 
upper limit drops to asi < 10~ 8 pb after imposing the constraints from colliders. In the 
region of the Higgs funnel where neutralinos are above 500 GeV, the spin-independent 
cross section drops to rather low values, asi > 7 x 10~ n pb, beyond the reach of future 
detectors. Basically the conclusions for the direct detection potential are similar to the 
ones reached for mt = 175 GeV, the best prospects are in the focus point region and for 
light neutralinos in the coannihilation region while the heavy neutralinos found in either 
the Higgs funnel or the coannihilation region will remain out of reach. 

To summarize the results in the mSUGRA model, the most promising models for 
both spin independent and spin dependent direct detection are models with a mixed 
bino/Higgsino LSP. For a top quark mass above 175 GeV, the former is expected only at 
large tan (3. Models with almost a pure bino LSP predict lower cross sections but should 
still be within reach of the large scale detectors provided neutralinos do not exceed roughly 
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Figure 15: Spin independent (top) and spin dependent (bottom) cross sections for scat- 
tering of neutralinos on proton as fonction of m^o for m t = 175GeV, tan (3 = 10 (left) 
and tan/3 = 50 (right). Black dots denoted models for which Qh 2 < .129. The line in the 
bottom figures is the expected limit from PICASSO. 
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Figure 16: Spin independent cross sections for scattering of neutralinos on proton as 
fonction of m^o for m t = 170 GeV, tan/3 = 10 (left) and tan/3 = 50 (right). Black dots 
denoted models for which Vth 2 < .129 



5.2 Non- universal models 

We next discuss non-universal models, considering the two cases that differ significantly 
from mSUGRA as far as the relic density is concerned: the case with a LSP with a higher 
Higgsino content, M3 < M2 as well as the case with a wino LSP, Mi > M%. Since the 
LSP coupling to the Z and/or Higgs is a crucial parameter in evaluating the neutralino 
proton cross section one expects significant differences with the mSUGRA model. 

5.2.1 M3 = M^/2 

First consider the model r^ = 2 and mt = 175GeV. We have performed a scan over the 
parameters M' 2 < 4TeV and Mo < 5TeV. As we have discussed before, both the squarks 
and the heavy Higgs tend to be lighter than in universal models. Furthermore light neu- 
tralinos have a higher Higgsino component since ji is typically smaller. Therefore one 
expects an enhanced scalar cross section. In FigElwe indeed see that in the coannihila- 
tion region the lower bound on the scalar cross section is around 5 x 10~ 9 pb and that it 
can exceed a si = 10 _8 pb for tan/3 = 10 for light squarks and neutralinos. The coanni- 
hilation/heavy Higgs region correspond to the inverted S-shape region in Fig. ITTk . The 
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enhancement of the cross section relative to the universal model is especially noticeable 
for heavy neutralinos. This is solely due to the enhanced coupling of the neutralino to 
the light Higgs. We had found two focus point regions in Fig. [TT] In Fig. [T7[ the first 
region corresponds to the sets of points around m^o pa 200 GeV for which asi ~ 10~ 8 , 
the second region to m^o pa 800 — 1200 GeV. Towards the upper range of this region 
the cross section drops to asi ~ 4 x 10~ 9 . As in the focus point region of the univer- 
sal model, reasonable cross sections are found despite the very heavy squarks because of 
the mixed bino/Higgsino nature of the LSP. There is finally one region where one finds 
a si ~ 10~ 9 pb, there neutralinos are light and annihilate through the light Higgs. We had 
already found a similar region in mSUGRA with roughly the same cross section. In this 
model, for tan/5 = 10, detectors such as Edelweiss probe part of the focus point region 
and almost all regions are within reach of future detectors like Genius. 

For tan/3 = 50, most of the parameter space is allowed by the relic density constraint. 
The predictions for asi spans a wide range of values. The largest cross sections are found 
in the low M — Mi/ 2 plane, however once the limit on the Higgs mass and on b — > 57 are 
taken into account one finds (J si < 8 x 10 _8 pb. Many models have asi < 10~ 9 pb, the 
absolute lower bound is 3 x 10~ 10 pb, an enhancement as compared to universal models. 
This is again directly related to the value of /i which entails a LSP with a higher Higgsino 
content than in mSUGRA. Some scenarios can be probed by EdelweissII but only if 
m^o < 400 GeV. To probe a large fraction of the parameter space one has to wait for large 
detectors like Xenon, even then it will not be possible to cover the full parameter space. 




Figure 17: Spin independent (top) cross sections for scattering of neutralinos on proton 
as fonction ofm^o for tan/3 = 10 and tanft = 50 in models where M' z = M' 2 /2. 
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5.2.2 M[ = 1.8M£ 



When M[ = I.8M2 and the LSP is wino-like, the region where the WMAP upper bound 
is satisfied includes all models with m^o < 600GeV. There the spin independent cross 
section ranges from agi ~ 10~ 9 — 5 x 10~ 8 pb for tan/5 = 10. For a given neutralino 
mass, the larger values for the cross sections, as usual, correspond to the lighter squarks. 
For heavier neutralinos, the WMAP allowed region is one of gaugino coannnihilation and 
heavy Higgs annihilation. Although coannihilation helps reduce the relic density it does 
not impact directly on the neutralino proton cross section, one finds small cross sections 
which further decrease as one increases m^o . Note that for a fixed value of the neutralino 
mass in the coannihilation/Higgs funnel region, the points that are allowed by WMAP 
have the lowest direct detection cross section since they correspond to larger values of M . 
For tan/3 = 10, detectors like EdelweissII can probe some of the models but even with 
Xenon, models with heavy neutralinos cannot be probed. For tan/5 = 50, the prospects 
for detection by large scale detectors are better, covering all the allowed parameter space. 
The main reason is that models with heavy neutralinos can only be found in the focus 
point region where the Higgs and Z exchange contribution is important. 
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Figure 18: Spin independent cross sections for scattering of neutralinos on proton as 



fonction of m^o for tan/3 = 10 and tan/3 = 50 in models where M[ = I.8M2. 



5.3 Summary 

In summary in the focus point region of mSUGRA the prospects are good for direct 
detection of dark matter by the next generation of detectors such as EdelweissII. In the 
stau coannihilation region the situation is more difficult and one will have to wait for large 
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scale detectors such as Xenon. Even then only a fraction of the parameter space can be 
probed. The Higgs funnel region is extremely dificult even for the most ambitious project 
as this region is also characterized by heavy neutralinos and rather heavy squarks, one 
has to face the problem of reduced sensitivity of detectors even when the cross section 
is not suppressed significantly. In models with non-universal gaugino masses M 3 < M 2 , 
prospects are better since one finds larger values for the neutralino couplings to the Z 
and the light Higgs. The lighter squarks also tend to increase the cross sections for a 
given neutralino mass. In models where M\ > M2 also the prospects for direct detection 
are more promising except in the coannihilation/Higgs funnel region at low tan/3. At 
large tan (3 the full parameter space can be basically covered since the models with heavy 
neutralinos and squarks that are hard to probe are excluded. 

6 Conclusion 

We have presented constraints on both mSUGRA and non-universal gaugino masses mod- 
els emphasizing the role of the latest measurements of WMAP on the relic density of dark 
matter. Our calculations are based on micrOMEGAs coupled to SOFTSUSY for the evalu- 
ation of the spectrum. 

We found, in agreement with previous analyses, that the new results of WMAP have 
almost excluded the bulk region at low M — M±/ 2 in mSUGRA and that it imposes very 
specific relation among sparticle masses and parameters: 

• Neutralino nearly degenerate with a slepton, 

• Neutralino mass near half the mass of a Higgs, 

• Neutralino with a significant Higgsino component. 

The latter is found in the focus point region although the location of the focus point 
region depends strongly on the value of the top quark mass and to a lesser extent on the 
choice of A . This comes as no surprise since the solution of the RG equation for fj, are 
very sensitive to the top Yukawa coupling. Furthermore the different implementations 
of the higher loop corrections in the RGE means that there are still differences between 
the predictions of different codes. In our case we found that a top quark mass above 
179 GeV pushes the focus point above M = 5 TeV. One is then basically left with 
the coannihilation region as well as a Higgs annihilation region. The latter includes a 
possibility for annihilation near a light Higgs resonance but annihilation through a heavy 
Higgs resonance is more likely to occur at large tan/3. Conversely if the top quark mass 
lies towards the lower end of the allowed range, a large focus point region is cosmologically 
allowed even at low tan (3. On the other hand the constraint from the light Higgs mass 
is more severe especially at low values of tan (3. Note that regions where the stop is the 
NLSP are usually ruled out by direct constraints, still the possibility exists in models with 
a large mixing in the stop sector. 

These conclusions result, in some sense, from the relation between the parameters 
of the MSSM obtained in the mSUGRA model. In order to assess how general the 
statements about the importance of the relic density constraints are, we have moved 



38 



away from mSUGRA by introducing some non-universality in the gaugino sector. By 
doing so one can deviate from mSUGRA predictions for the parameters that are critical 
in the evaluation of the relic density in particular the value of fi. This parameter in 
relation to Mi,M% determines the gaugino/Higgsino content of the LSP and as a result 
directly the annihilation cross sections into fermions and gauge bosons. The impact of 
this parameter is already visible in mSUGRA where as soon as one could get away from 
the pure bino case, one finds (in the focus point region) much more efficient annihilation of 
neutralinos. In non-universal MSSM we have, in addition to a Higgsino LSP, models with 
a wino LSP. In these models the upper bound on the relic density can easily be satisfied, 
because of the large cross section for the annihilation of neutralinos into W pairs. These 
models are the ones where M[ > M' 2 at the GUT scale. Second, by changing the relation 
between the heavy Higgs mass and the neutralino mass we have found that annihilation 
through a heavy Higgs exchange can take place even at low tan (3. These conditions can 
be realized in non-universal gaugino mass models, by changing the mass relation between 
M3 and M' 2 /M[. Indeed Ma is sensitive to the value of M' z at the high scale. Finally in 
all models we found a region where the WMAP upper limit could be satisfied because 
of the contribution of slepton coannihilation channels as long as the neutralino and stau 
were not too heavy. In summary in non-universal gaugino mass models, compatibility 
with the WMAP measurement imposes similar conditions on the supersymmetric masses 
and parameters than in mSUGRA. The only new possibility is for a wino LSP in addition 
to the Higgsino LSP. We stress that SUGRA models with non-universal gaugino masses 
give a good picture of the dark matter constraints in the general MSSM while keeping 
with a reasonable number of free parameters. Although we have mainly studied only two 
non-universal gaugino mass models, it was sufficient to illustrate the main mechanism at 
work in the relic density. 

We have also discussed the potential for probing SUGRA models with the different 
experiments for direct detection of the LSP. While few models lead to prediction for the 
direct annihilation cross section within reach of the current experiments, a large fraction of 
the parameter space will be probed by ton-scale future detectors. Typically models with a 
Higgsino LSP predict the largest cross sections while coannihilation models feature much 
lower cross sections, especially when the sfermion masses increase. In fact the Higgsino 
component of the LSP is also what makes efficient annihilation of neutralino pairs and 
leads to acceptable values for the relic density. Finally when the neutralino annihilates via 
a heavy Higgs, few models can lead to a signal for direct detection especially when they 
feature a heavy LSP. We keep a more detailed comparative study of the direct detection 
and collider searches for a future publication. 
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